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Abstract. We study the wave equation in a bounded domain 
or on a compact Riemannian manifold with boundary. Assume 
that we are given the hyperbolic Neumann-to-Dirichlet map on 
the boundary corresponding to physical boundary measurements. 
We consider how to focus waves, that is, how to find Neumann 
boundary values so that at a given time the corresponding wave 
converges to a delta distribution S y while the time derivative of 
the wave converges to zero. Such boundary value are generated by 
an iterative sequence of measurements. In each iteration step we 
apply time reversal and other simple operators to measured data 
and compute boundary values for the next iteration step. The key 
feature of the algorithm is that it does not require knowledge of the 
coefficients in the wave equation, that is, the material parameters 
inside the media. However, we assume that the point y where the 
wave focuses is known in travel time coordinates. 
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1. Introduction 
Let us consider the wave equation in a bounded domain M 

{u u (x, t) + Au(x, t) = 0, in M x 
u\ t =o = 0, u t \ t =o = 0, 
dv u \dMxR+ — f: 



- j 



where A is a 2nd order elliptic partial differential operator. 

In this paper we show how to construct Neumann boundary values 
/ such that at time T, the wave (u(T),u t (T)) is arbitrarily close to 
(cS y , 0) , where 5 y is the Dirac delta distribution at a chosen point 
y G M . We call such waves focusing waves. To find such boundary 
values, we only assume that we can make physical measurements from 
the boundary of M . For given Neumann boundary values we can 
measure the Dirichlet boundary values of the wave. A focusing wave 
can then be generated by an iterative sequence of measurements. In 
each iteration step we apply time reversal and other simple operators 
to measured data and compute boundary values for the next iteration 
step. 
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The iteration algorithm in this paper is closely related to time rever- 
sal methods. Let us therefore shortly discuss the underlying idea and 
the usually used approximations behind these methods. As a simple 
example, let us consider a domain M in M 3 , and suppose that we can 
measure waves and generate sources on the boundary of M . Let us 
first assume that there is a theoretical point source at y G M , and 
we measure the wave and/or its normal derivative at the boundary of 
M . Assume further that we record this signal, reverse it in time, and 
re-emitted into M , see [19]. Then one can show (assuming certain 
approximations hold, see [TH, EH EE]), that the re-emitted wave will 
travel like the original wave, but as if time were running backwards. 
This causes the re-emitted wave to focus near y . 

This principle can also be used for imaging. To find a small scatterer 
D in a relatively homogeneous domain M , one sends a wave into 
M . If the scatterer is small and the single scattering approximation 
is justified, the scattered wave corresponds to a wave produced by a 
point source at D . If we record this scattered signal at the boundary, 
reversed it in time, and re-emit it into the domain, it will focus at the 
scatterer. Furthermore, this focusing has been observed to be quite 
stable under perturbations of the medium. Thus, if the re-emitted 
wave is simulated (by computational means) in homogeneous media, 
it will focus at the location of D. In this way a small scatterer can 
be found using relatively simple computational methods. The above 
basic idea has been refined in various ways. If the target area contains 
multiple scatterers, an iteration scheme can be used to focus the wave 
on any of the scatterers [33] . 

Besides imaging, time reversal can be used to focus a wave onto 
a scatterer, say, inside the human body. One application of this is 
litotripsy, where one breaks down a kidney or bladder stone using a 
focusing ultrasonic wave. Another application is hyperthermia, where 
a cancer is destroyed by an excessive heat dose generated by a focusing 
wave. Let us point out that for the wave equation, there are various 
methods to estimate material parameters in travel time coordinates 
from boundary measurements. These methods are, however, quite un- 
stable [21 EH]- Therefore they might not be suitable for hyperthermia, 
where safety is crucial. An important question is therefore how to focus 
waves in unknown media. 

For reviews and extensions on time reversal, see seminal papers of M. 
Fink, [dEElED]- Time reversal methods have been intensively stud- 
ied in random heterogeneous media where the statistic of the random 
media is known, see e.g. [51 EH 12 [121 E]- For time reversal in chaotic 
cavities, see [39]. For related analysis on time reversal methods, see 

also [aniiisiEniEi]. 

Let us describe the key features of the algorithm in this paper. First, 
to focus a wave we do not require knowledge of the material parameters 
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inside the media. We only assume that the coordinates of the point 
are known in travel time coordinates. This means that focusing can 
be done in the same coordinates in which imaging is done. Thus, as 
the algorithm for focusing does not require media parameters obtained 
from imaging, errors in imaging do not accumulate into errors in fo- 
cusing. Second, the algorithm can focus waves near an area having 
no scatterers. Third, the algorithm is computationally cheap. In a 
sense, all computations are done in the media; there is no need to solve 
the wave equation, cf [23]. We will assume that the medium is linear, 
non-dispersive, non-dissipative, and frequency-independent. However, 
we do not need any other approximations like the single scattering 
approximations to prove that the algorithm works. 

The limitations of the present algorithm is that we assume selfad- 
jointness of operator A and that time T is large enough. Moreover, 
the point y where the wave focuses need to be specified in travel time 
coordinates unless operator A is known. 

The present work is a continuation of [UJ where a similar iterative 
scheme was introduced, for which u{T) focuses to a delta-distribution, 
but the time derivative u t (T) is uncontrolled. The present work can 
also be seen as also a generalization of so-called retrofocusing in control 
theory, where the aim is to produce boundary sources giving the same 
final state as a boundary sources sent before in the medium, see [241IBT]. 
The methodology in this paper arises from boundary control methods 
used to study inverse problems in hyperbolic equations [HI El QUI EH 



The outline of this work is as follows. In Section El we introduce 
notation and review some relevant results from control theory. We also 
define the boundary operators that are needed in the iteration scheme. 
In Section [3] we describe the main results (Theorems [5] and [7]) and 
outline their proofs, and in Section [H we prove these results. 



We assume that M C M m (m > 1) is the closure of an open 
C°° -smooth bounded set with non-empty smooth boundary dM or 
a C°° -smooth compact manifold with boundary. Furthermore, we 
assume that M is equipped with a C°° -smooth Riemannian metric 
g = ^2jk 9o k ® dx k . Elements of the inverse matrix of gij are 
denoted by g % i . Let dV g be the smooth measure 



where \g\ = det([gjk\) ■ Then L 2 (M) is defined by the inner product 



EZlEElES]. 



2. Definitions 
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g(x)\ 1 ^ 2 dx 1 A ■ ■ ■ A dx m , 
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where dV = fidV g and fi G C°°(M) is a fixed strictly positive function 
on M . 

In wave equation (Dp), we assume A represents the most general 
formally selfadjoint elliptic partial differential operator with respect to 
the above inner product [27]. In local coordinates, A has the form 

m 1 d ( dv\ 

*• = -E Mx)|gM |i/ 2 ^j {K'M'WWtoi ) + 9(*)«. 

where g is a smooth function g: M — > K. For example, if \i = 1 
and g = then .4. reduces to the Riemannian Laplace operator. Let 
us point out that A represents media that is linear, non-dissipative, 
non-dispersive, and frequency-independent. 
On the boundary, operator d v is defined by 

m d 
i=i 

where v(x) = (v 1 , z/ 2 , . . . , z/ m ) is the unit interior normal vector satis- 
fying YlTk=i9jk l/: 'v k = 1 . To integrate functions on dM we use the 
measure dS on dM induced by dV g . If B C dM x R + , we define 

L 2 (B) = {/ G L 2 (<9M x R+) : supp(/) C B} 

identifying functions and their zero continuations. 

With these assumptions, the wave equation has a solution whenever 
/ G L 2 (dM x M+) , and we denote this solution by . The map 
/ i — ► is linear over R , and = m 9 *-^ when /, 9 t / G L 2 (dM, R + ) . 

Let <i(x, y) be the geodesic distance corresponding to g . The metric 
d is also called the travel time metric because it describes how solutions 
to equation (OQ) propagate. By the finite velocity of wave propagation, 
(see [22]) we have that if V C dM is open, and / G L 2 (T x E + ) , then 
at time t > , solution it-' is supported in the domain of influence 

M(T,t) = {xeM : d(z,T) < 0- 

The diameter of M is defined as 

diam(M) = max{d(x, y) : x, y G M}. 

The characteristic function of a set 5 is denoted by xs ■ 

2.1. Controllability for wave equation. The seminal result imply- 
ing controllability is Tataru's unique continuation result [35l 137] . 

Proposition 1 (Tataru). Let u G Hl oc (M x E + ) 6e a solution of the 
wave equation 

u tt {x, t) + ^4n(x, t) = 0. 

Assume that 

w|rx(o,2r) = 0, d u u\r x (o,2r) = 0, 
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where T C dM is an non-empty open set and r > . Then 
u(x, t) = 0, d t u(x, t) = forxE M(T, r). 

Using Tataru's unique continuation result, one can prove the follow- 
ing controllability results (The proof is postponed to section H]): 

Proposition 2 (Approximate global controllability). If T > diam(M) , 
then the linear subspace 

{V(T),uf(T)): /GC °°(9MxR + )} 

is dense in if^M) x L 2 (M) . 

This result yields the following controllability result, see e.g. [27] 
and references therein. 

Proposition 3 (Approximative local controllability). let r > 0, let 

Fx, . . . , T j C dM be open non-empty sets, and let < Sk < r , k — 
1, . . . , J . Suppose 

j 

B = \jY j x (r- Sj ,r), 
i=i 

and P is multiplication by the characteristic function xb , 
(2) P: L 2 (dM x (0,r)) -> L 2 (dM x (0,r)), 

f(x,t) ^ Xjs(x,t) f(x,t). 

Then the linear subspace 

{u Ph {r) : h E L 2 (dM, (0,r))} 

is dense in L 2 {N) , where N = Uj=i M(Tj, Sj) . 

2.2. Operators for boundary sources. In this section we introduce 
operators for manipulating boundary sources. These will be needed 
both in the proof of the main result and in the iteration scheme. 

For initial boundary value problem (pQ) we define the non-stationary 
Neumann-to-Dirichlet map (or response operator) A by setting 

Af = u f \ dMxR+ , feL 2 (dMxR + ). 

In other words, we solve the wave equation (CD) for a boundary source 
/ , and measure boundary values for the solution when t > . In 
this work, we only need the finite time Neumann-to-Dirichlet map, 

-A-2t/ = u ^\dMx(0,2T), 

where T > . By |36| the map 

A 2T : L 2 (dM x (0,2T)) -> H 1/3 (dM x (0,2T)) 

is bounded, where H s (dM x (0,2T)) is the Sobolev space on dM x 
(0,2T). 
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For / E L 2 (dM, (0,2T)) , let 

Rvrffat) = f(x,2T-t), 

J2rf(x,t) = / Jzr(s,t)f(x,s)ds, 

J[0,2T] 

where J 2 t(s, t) — t) and 

L = {(s,t) E M+ x R + : t + s < 2T, s > t}. 

We call R2T the time reversal map, and J2T the time filter map 1 1 1 J . 
On L 2 (<9M x [0,2T]) with the measure dS{x)dt , the adjoint of A 2T 
is dU, 

Agy = R2T^-2tR2T- 

For / e L 2 (<9M, [0,2T]) , let 

Qirf = / g(t,s)f(x,s)ds, 
Jo 

be the time filter operator, where g: (0,2T) 2 — > JR., 

, . 1 f(e _t + e i )(e 4T e~ s + e s ), i < s, 

!/(M) ~ 



2(e 4T - 1) \(e- s + e s )(e 4T e-* + e<), t > s, 

is the Green's function for the problem 

l-d 2 )g(t,s) = 5(t-s), 
3t0|t=o = O, %|*= 2 t = 0, se(0,2T). 

Next we consider A 2T , R2T, J2T, Q2T as operators such that 

A 2T , i? 2T , J 2T , Q 2T : L 2 (dM x [0,2T]) -> L 2 (3M x [0,2T]). 

Below, we often denote _R 2 t, J2T, and Q2T by i?, J, and Q . For 
f,h E L 2 {dM x [0,2T]) the Blagovestchenskii identity states that 

f u f {T)u h {T)dV = [ (Kf)(x,t)h(x,t)dS(x)dt, 

JM JdMx[0,2T] 

where K: L 2 L 2 is the bounded operator 

K = K2T — -R 2 tA 2 T-R2T<^2T — J 2 tA 2 t- 

For a proof, see e.g. [11]. The importance of this identity is that 
it shows that the inner product of solutions u^{T) and u h {T) can 
be calculated from the boundary. Namely, on the right hand side of 
the Blagovestchenskii identity, dS is the Riemannian surface volume 
on dM , and K is defined in terms of the Neumann-to-Dirichlet map 
A 2 t and simple operators on boundary values like time reversal. The 
intrinsic Riemannian surface volume dS on dM is determined by A 2 t . 
Namely, by Tataru's unique continuation principle, the Schwartz kernel 
of A 2T is supported in 

E — {(x, t, x\ t') E {dM x [0, 2T]) 2 : t - t! > d{x, x')}, 
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and the boundary dE is in the support. The set dE determines the 
distances of points z, z' G dM with respect to the intrinsic metric of 
the boundary (dM, gahi) ■ 

3. Iterations and main results 

3.1. Cutoff of wave. In this section we describe Theorem [5] which 
can be seen as a lemma used in the proof of Theorem [7J 
Let X be the Banach space 

X = L 2 (dM x [0, 2T]) xY, Y = H 1 ^, 2T); L 2 (dM)), 

such that the inner product on X is 

r hi\ fh 2 



CL\ I \a 2 



= (hi, h 2 ) L * + (ai, a 2 ) L 2 + (d t cti, d t a 2 } L 2 - 

x 

Definition 4. Let T > 2 diam(M) , and let 

j 

B = |Jr j x(r-^ ) 

i=i 

where Ti,...,Tj C dM are open non-empty sets, and < Sj. < T, 
k = 1, . . . , J . Let P = Xb be the multiplication with the characteristic 
function of B defined as in equation (J2J), and let L : X — > X be the 
operator 

'I 0\ (2PKP -PK \ 
Qj {-KP K-dtKdt) ■ 

Let a G (0, 1) . and let uj > be such that 2(1 + ||£||x) < oj , and let 

uj uj 

If f G L 2 (dM x E + ) 6e a boundary source, we define a sequence 
' 61, n = 1,2,... % 



(3) 




(4) 



Theorem 5 (Cutoff of wave). Let ai(a), a 2 (a), . . . , be as in Definition 
[|| Then the sequence converges in Y , 

lim a n (a) = a(a), 



and function a(a) G F on the right hand side satisfies 

/„««»(T)\ _ ( x „uf(T)- 
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where both limits are in L 2 (M) and N is the domain of influence 



J 



n= |jM(r fc , Sfc )- 



fc=l 



Note that here lo may depend on a . For instance, we can choose 
to — 1/a . 

Let us emphasize that the novelty of this theorem is the explicit itera- 
tion scheme for a(a) depending only on boundary measurements. The 
scheme depends on operators J,Q,P , and K that can be calculated 
from the boundary of M . The first three are simple operators like 
integration and restriction. Operator K = Rit^itRitJit — Jit^-it 
involve time reversal R<it , time filtering Jit and two evaluations of 
the Neumann-to-Dirichlet map A 2 t which corresponds to two physical 
measurements. Hence, the first order approximation of a(a) requires 
2 physical measurements. After that, each additional term requires 10 
additional measurements. Thus, for finite approximation of a(a) we 
only need finitely many evaluations of the Neumann-to-Dirichlet map. 

The full proof of Theorem [5] is given in Section 14. 1L Let us here 
outline the main ideas. For a G (0,1), boundary sources h(a),a(a) 
are defined as the minimum of the functional 



In the sequel, when there is no danger of misunderstanding, we denote 
the L 2 -norms in the spaces L 2 (M) , L 2 (dM x [0,T]) etc. just by 
|| ■ || . In LemmaEl we use convexity to prove that for each a , there is a 
unique minimum h(a), a(a) , and by studying the Frechet derivative of 
T we find a linear equation (equation (jHJ)) for this minimum. In Lemma 
M we show that iteration scheme (p!J converges to a von Neumann sum 
that represents the solution to equation (jHJ). That minimizer a(a) 
satisfies the sought limit is proven in Lemma [101 The key step in the 
proof is to use the approximative controllability results from Section 
12.11 to show that the first terms in T can be arbitrarily close to ||(1 — 
XaO^(T)|| 2 and the next two terms can be made arbitrarily small. 

3.2. Focusing of wave. To understand how one can focus waves us- 
ing Theorem [5] suppose we have sets B C B C dM x [0, T] (defined in 
terms of Tj and Sj as in Definition BJ . Then Theorem [5] implies that 



T{h, a, a) 
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there are boundary sources a(a) and a(a) such that 



lim 



u 



° {a HT)\ (xnuHT) 



where the domains of influences satisfy N C N at time T . As the 
solution operator / i— > is linear and commutes with <9 4 , solution 

6(a) = a(a) — a (a) satisfies 



lim 



« 6(a) cn\ / XN\N uf (r) 

ut\T)) \ 



That is, in the limit, the solution corresponding to b is supported in 
N \ N . In the proof we construct P and P such that N \ N is a 
family of sets that shrink onto a chosen point x . By further scaling b 
with a suitable constant depending on the volume of N\N , we obtain 
the delta distribution. 

To formulate the result, let us introduce some notation. By 7^ we 
mean a geodesic in (M, g) parametrized by the arclength such that 
7*.f(0) = x, 7^(0) = C and = 1 . Let 1/ = z/(^) , z e be 
the interior unit normal vector to dM . Then there is a critical value 
t(z) > , such that for t < t(z) the geodesic 7« )I /([0, i]) is the unique 
shortest geodesic from its endpoint J z , v (t) to dM , and for t > t(z) 
it is no longer a shortest geodesic. We will not consider the degenerate 
case t = t(z) . 

Definition 6. Let T > 2 diam(M) , let x — 7z^(T) , where z G dM , 
and < T < T . Let Tj C dM for j = 1, 2, . . . be open sets around 
z , such that Tj D and HJli = {^} ■ 

Suppose f G C , ^°(9MxIR + ) . Let a n (a,e) G F be functions obtained 
from the iteration in Definition^ when B is the set 



B(e) = dM x (t- (f -e),T 

a G (0, 1) , and e > is sufficiently small. Similarly, let a n (a,j,e) G 
Y be functions obtained from the iteration in Definition^ when B is 
the set 

B(j, e) = {dM x (T - (f - e), r)) U (r, x (t — (T + e),T 

a G (0, 1) . j — 1, 2, . . . and e > is sufficiently small. 
Under these assumptions, let 

bn(a,j,e) = (a n (a,j,e) - a n (a,e)) G Y. 

Theorem [7] is the main result of this paper. 
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Theorem 7 (Focusing wave). Let z G dM , T, b n (a,j,e) for n = 
1,2,... be as in Definition® Then functions bj converge in Y , 

lim b n (a, j,e) = b(a,j,e). 

n— >oo 

Moreover, if T < r(z) then functions b(a,j,e) G Y satisfy 

(6) (^(T)) = C ® uf ( T >*) (o) ' 

where the inner two limits are in L 2 (M) and the outer limit is in 
3t'(M) , and C(x) is defined in equation (fl~3l) . 
If T > t(z) , the limit (0|) is zero. 

Let us make three comments about this theorem. First, the condition 
T < t(z) means that z G dM is the closest boundary point to x 
and d(x, z) = d(x, dM) . If a point x has a unique closest boundary 
point, we say that it is admissible. For example, on the closed disc, all 
points are admissible except the center. For a general (M, g) the set 
of points that are not admissible has measure zero. Second, we assume 
that / G C^{dM x R+) . Hence u f G C°°(M x M+) (see [32]), and 
uf(x,T) exists pointwise. Third, a function v G L 2 (M) is interpreted 
as a distribution v G 3>'{M) by the formula 

(v,(j>) = [ v<f)dV, G @(M). 

The delta distribution at y G M is defined by (5 y , 0) = <fi(y) for 
(j) G @{M) . 

4. Proofs 

We start with the proof of Proposition [2l The proof is a relatively 
direct consequence of Tataru's unique continuation theorem and can be 
found e.g. in the case of Dirichlet boundary conditions in [3H Lemma 
2.1]. 

Proof of Proposition [1 Assume that a pair 

(ip, -0) G (H\M) x L 2 (M))' = Hq\M) x L 2 (M) 
satisfy the duality 

(u f (T),i/j) iHl{MlH -i {M)) + (u{ (T), -4>) L 2 (M) = 

for all / G C °°(9M x (0,T)) . Note that H\M) is the domain of the 
square root of the operator A + cl when c is large enough, denoted 
by V(A 1/2 ) and Hq\M) is the dual H\M) = V{A~ l/2 ) . Let 

e tt + Ae = inMx(0,T), 

d u e\dMx(o,T) = 0, e\ t=T = <fi, e t \ t=T = ip. 
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By [32], e G C\[0, T], L\M)) n C([0, T], ^(M)) and eUx(o,T) G 
H~ 2 ' 5 ~ e (dM x (0, T) , e > . Thus we have in sense of distributions 

= / [w 7 (e« + -4e) - (u£ + At^e] dt 

JMx(0,T) 

u{ {T)(j)-u s {T)i))dV + f fe dS x dt 

M JMx(0,T) 

— I f e dS x dt 

JMx(0,T) 

for all / G C™(dM x (0,T)) . This yields that 

e|dA/x(o,T) = d v e\QMx(o,T) = 0. 

To apply unique continuation for e G C([Q,T], Hq 1 (M)) , let e > 0, 
and let 77 G Co°(R) be a function supported on (—1,1) C R whose 
integral over R is one. Then 

/" £ - t' 

e e (a;,t) = / e{x,t')r)( )dt' 

satisfies 

(d 2 + A))e e = in M x (e, T - e), 9„e e | aif x( E ,T- £ ) = 

and e e G C°°((e, T— e), Hq 1 (M)) . By representing e t in terms of eigen- 
functions of A, we see that e e G C°°((e,T - e),V{A°°)) C C°°(M x 
(e, T — e)) . Using Tataru's unique continuation theorem [35] we see 
that if < e < (T — diam (M))/2 then e £ (T/2) = d t e t (T/2) = 0. 
Hence e e = identically on M x [0, T] . When e — > , we see that 
also e vanishes identically and thus (ft = ip = . □ 

4.1. Proof of Theorem [5l On X we will study the minimization 
problem 

(7) min T(h , a, a), 

(h,a)GX 

where a G (0, 1) and T is defined in equation ([S]). By [32], the map 
i — ► is continuous L 2 (dMx [0,T]) -> C([0, T]; iJ 5 / 6_e (M)) , e > 0. 
Thus (h, a) h- > .F(/i, a, a) is continuous map X — > R . 

Lemma 8. For any a G (0, 1) minimization problem ([7)) has a unique 
minimizer (h, a) G X . This minimizer is the unique solution to 



(8) (a + L) 



& - n 



where L is defined in equation ([3]) . Furthermore, L : X — > X is non- 
negative, bounded, and self adjoint. 
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Proof. We have 

F(h,a,a) = (f - Ph, K(f - Ph)) + (Ph - a, K(Ph - a)) 
+ (d t a, Kd t a) 

+ a((h, h) + (a, a) + (d t a, d t a)), 

Here K and P are selfadjoint in L 2 (M) . The dual F* of the Hilbert 
space Y = H l {[0, 2T}; L 2 {M)) can be identified with H~ l {[0, 2T};L 2 (M)) 
As Q the inverse of 1 — d 2 , it can be considered as an operator 
Q : Y* — > Y . Thus using d t Kd t : F — > F* , we can write 

(d t a,Kd t a) L 2 {M) = -(d t Kd t a,a) Y *,Y = -(Qd t Kd t a,a) Y - 

Thus, using the inner product on X = L 2 (M) x F , we can rewrite JF 

as 



F{h,a,a) = (f,Kf) + 2 



h\ (1 0\ f-PKf 
a J ' 10 QJ\ 



A 



As Q: F* -> F and <9 t ifd t : F -> F* are bounded, L: X -> A is 
bounded. A direct calculation shows that L is self-adjoint. Setting 
/ = and a = in equation Q shows that L is non-negative. 

Let us next observe that T is strictly convex, so for the first claim 
it suffices to prove existence; by convexity, a local minimum is a global 
minimum, and by strict convexity, there is only one global minimum 
(see [311 Theorem 38. C]). To prove existence, we recall that by [T], 
x G X is a local minimum of T : X — »• K provided that: 

(i) the first two Frechet derivatives at x , D X T{£) and D 2 X T{^, rf) , 
exist and are continuous, 

(ii) D x F = 0, 

(iii) u t— > D^J 7 ^, ■) is a linear isomorphism X —> X* , and 

(iv) Z&F(£, > for all f ^ . 

The sought Frechet derivatives of T : X — > R are 

«) ("T / ) +(a+i) C)' 5 

D 2 J~h, a (£, Tj) = 2(Z,(a + L)r ] ), £, V G A. 
It remains to prove that a + L is invertible, but if (a + L)(£) = , 
then by non-negativity, < (L£,x)x = — a \\C\\x — 5 an d £ = • ^ 

Lemma 9. Iteration scheme (J4]) converges to the unique solution to 
equation (jHJ). 

Proof. Using S 1 and u defined in definition [4J we may rewrite equation 
(ED as 



(I-S) 



a/ 
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For a self-adjoint operator B , the norm satisfies ||B|| = supjj (x, Bx) | : 
||x|| = 1} . Hence a + (Lx, x) < u/2 , and by non-negativity, 

\\S\\<1-™<1. 

UJ 

We may therefore iteratively solve h, a by a convergent von Neumann 
sum. □ 

Lemma 10. Minimizers h(a),a(a) G X for ([7]) satisfy 

supp h(a) C B, 

a(a) G range Q, 
/«■*■> (T)\ (xNuHTf 

where all limits are in L 2 (M) . 

Proof. The first two claims follow by writing out ((HJ). For the other 
results, let us define Z : X — > E by 

Z{h,a) = \\\xNU f {T)-u p \T)\\ 2 + \\\u a {T)- XN uf{T)\\ 2 

+ IKCOII 2 . 

To prove the last claim we show that for any e > there exists an 
a(e) G (0,1) such that Z(h(a), a(a)) < 4e when a G (0,a(e)). Let 
us note that 

u p f(T) = XNU Pf {T), f G L\dM, [0, T}). 

Hence, for any (h, a) G X , 

F(Ph,a,a) = \\(l-XN)u f (T)\\ 2 + \\ XN uf(T)-u Ph (T)\\ 2 

+ \\u Ph {T)-u a {T)\\ 2 + \\u a t {T)\\ 2 

+ a(\\h\\ 2 + \\a\\ 2 + \\d t a\\ 2 ). 

It follows that for any (h, a) G X and a G (0, 1) , 

Z(PM) < F(Ph,a,a)- - XnW (T)\\ 2 . 

Here we have estimated the second term in Z using the triangle in- 
equality and the inequality (s + 1) 2 < 2(s 2 + 1 2 ) . Let us fix e G (0, 1) . 
By Proposition [3] there exists an h £ G L 2 (B) such that 

\\ X NU f (T)-u Ph *(T)\\ 2 <e, 

and by Proposition [2] there exists an a £ G H 1 ^) such that 

\\u a *{T)- X NU Ph *{T)f < e, 
|K(T)|| 2 < e. 

As /i e = P/i e we have 

F(h £ , a £ , a) = ||(1 - Xtv)^(T)|| 2 + 3e + a (||/i e || 2 + ||a £ || 2 + \\d t a £ \\ 2 ) , 
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and if a G (0, 01(e)) , where 



ale) 



l+IIMI 2 +lk|| 2 + ||<9 4 a £ || 2 ' 
then the minimizer h(a),a(a) of T satisfies 

Z(h(a),a(a)) < F(Ph(a), a(a), a) - ||(1 - X N)u f (T)\\ 2 

< F(h e ,a e ,a)-\\(l- XN )u f (T)\\ 2 

< As. 



□ 



4.2. Proof of Theorem [S Let us note that for any e > , then 
T k C B(z,e) for sufficiently large k. By Theorem [5], the following 
limits exist in Y , 

a(a,e) = lim a n (a,e), 

n— >oo 

a(a,j,e) = lim a n (a,j,e). 



<^(T)\ (xN { s)uf(T) 
f a > e \T)J V 

(aJ ' £) (T)y V 



N(e) = M(dM, T — e), 
N(j,e) = M(dM,f -e)UM(T j ,f + e). 
We define b(a,j,e) = lim n _ >00 6 n (a, j, e) , whence 

b(a,j,e) = e *~ (a(a,j,e) — a(a,e)) . 
Lemma 11. In L 2 (M) , 



and 






(10) 


lim 






\u 


(11) 


lim ( 


' u a{ 
a( 






where 







t ^) (T )X _ / XJ(£)M /(T) 



where 

J(e) = M(z, T + e)\ M(dM, f-e). 
Proof. Since ahh' is linear, it suffices to prove that pointwise 

j l™X M (r J ,f +£ )\M(9M,f- £ )( I ) = XJ(e)(x), xeM. 

This is clear for x G J(e) . If x ^ J(e) we claim that x ^ M(T,,-,T + £:) 
for large j . However, if d(x, z) > T + e , then 
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for large I . For y G Ti , d(x, y) > d(x, z) — d(y, z) > T + e , so 
d(x,Fi) > f + e, and x £ M(T h f + e) . □ 

The next Lemma show that J(e) are sets that shrink onto x in the 
case when T < t(z) . 

Lemma 12 (Properties of J(6)). For any e > , there is a 5 > 
such that 

(12) J(5)cB(x,s). 

Moreover, if T < t(z) , then {x} C J (6) for all 5 implying that 
n<5>o ^(^) = {^} ■ If T > t(z) then J(S) = for 5 small enough. 

Proof. Let us first prove (fT2l) . For a contradiction suppose that e > 
and X\,X2, ■ ■ ■ is a sequence such that 

XjeJ(l/j), Xj£B(x,e). 

As M is compact, we can move onto a subsequence and assume that 
Xj converges to an x G M \ B(x,e) . Now d(xj,z) < T + 1/j and 
d(xj, dM) > T — 1/j , and as d(x, dM) is continuous, 

d(x, z) < f, d(x, dM) > f. 

Thus f < d(x, dM) < d(x, z) < f . If f > r{z) , we have d(x, dM) < 
T and obtain a contradiction. Thus we can assume that T < r(z) . 
Then the above inequalities yield that d(x, z) = d(x, dM) = T . As M 
is compact, there is a geodesic from z to x that realizes d(x, z) . Then 
7] also realizes d(x, dM) , and r] must necessarily be normal to dM 
[H]. Thus rj — 7 , and contradiction. Thus (fT2l ) is proven. 

If T < r(F) , then clearly x G J(5) for all 5 > . On the other hand, 
if T > t(z) then there is G (9M such that d(x, z') < d(x, z) and 
we see that x G" J(S) for small <5 . Thus we have shown that J(5) = 
for 5 small enough. □ 

Proof of Theorem^ Consider first the case when T < t(z) . Then we 
observe that the following limit exists 

exists (see [H]). Here Yo\(A) = J A ld\^ when A C M . Let us also 
note that B(x,e/2) C J(e) so Vol(J(e)) > 0. Thus, as m^(T, •) is 
continuous, 



^oVol(J( £ ))y j( 



= (C@)u f (T,x)6s,<l>), (f>e®(M). 
The result follows by [ZD, Theorem 2.1.8]. 
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In the case when T > t(z) , J(e) = for e small enough, and thus 
the limits (fTUl) and (fTTl) are the same. Hence limit (jSJ) is zero. □ 
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